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Dark matter annihilation or de-excitation, decay of metastable species, or other new physics may 
inject energetic electrons and photons into the photon-baryon fluid during and after recombina- 
tion. As such particles cool, they partition their energy into a large number of efficiently ionizing 
electrons and photons, which in turn modify the ionization history. Recent work has provided a 
simple method for constraining arbitrary energy deposition histories using the cosmic microwave 
background (CMB); in this note, we present results describing the energy deposition histories for 
photons and electrons as a function of initial energy and injection redshift. With these results, 
the CMB bounds on any process injecting some arbitrary spectrum of electrons, positrons and/or 
photons with arbitrary redshift dependence can be immediately computed. 
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I. INTRODUCTION 

Precise measurements of the cosmic microwave back- 
ground (CMB) anisotropies (particularly by WMAP [1 
and the upcoming Planck Surveyor [2] ) have the potential 
to probe the detailed ionization history of the universe 
between recombination and reionization. In particular, 
dark matter (DM) annihilation or similar new physics can 
inject high-energy electrons and photons into the photon- 
baryon fluid, which give rise to additional ionization and 
heating as they cool, and consequently leave an imprint 
in the CMB [3H5] . 

High-energy electrons in the early universe dominantly 
lose their energy by upscattering CMB photons, and the 
cooling time for photons - either upscattcrcd or injected 
directly - can easily be comparable to a Hubble time 
[5HZ]. Consequently, the effect on the ionization history 
and gas temperature does not precisely track the energy 
injected; particles injected at a particular redshift may 
lead to an increase in the residual ionization fraction at 
a considerably later time. 

In this note we employ code developed to map 
an energy injection history into an energy deposi- 
tion history, and described in detail in 5 . Wc 
compute the deposition histories for a large grid of 
photon and electron injection energies and redshifts, 
and make the results publicly available online, at 
http : //nebel . rc . f as . harvard . edu/ epsilon; interpo- 
lation on this grid can be used to compute the energy 
deposition history for an arbitrary spectrum and varia- 
tion of the energy injection with redshift. 

We describe our numerical results in Sec. |Hj and dis- 
cuss their interpretation in Sec. |III| as a sample ap- 
plication, we demonstrate the effect of including DM 
halo formation at low redshifts on the energy-deposition 
history. In Sec. |IV[ we briefly review the principal- 
component method for estimating CMB constraints de- 
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veloped in |23j : in Sec. |vj we outline some applications, 
showing as examples how our results can be applied to 
set new and updated constraints on scenarios with late- 
decaying species, and models of DM annihilation where 
the annihilation switches on with some characteristic long 
timescale. In Sec. VI we present our conclusions; Ap- 



pendix [X] describes the online supplemental material. 

Throughout this work, we assume the cosmological pa- 
rameters from [5]: explicitly, ujb = 2.258 x 10~ 2 , lj c = 
0.1109, ^(£=0.002 Mpc" 1 ) = 2.43 x 10~ 9 , n s = 0.963, 
t = 0.088, H = 71.0 km/s/Mpc. 



II. DESCRIPTION OF NUMERICAL RESULTS 

The code developed for [5] takes as an input some in- 
jection of photons and electrons, with a specified redshift 
and energy dependence. Backreaction on the CMB pho- 
tons and gas is not included, as large modifications to 
the ionization history or CMB spectrum are ruled out by 
observational constraints, and consequently the problem 
is (to a good approximation) linear. We thus populate 
individual energy bins with electrons/positrons or pho- 
tons at a specific redshift and track the spectral evolution 
with redshift. Our 40 energy bins are log-spaced between 
1 keV and 10 TeV, in photon energy and electron kinetic 
energy. We employ 65 log-spaced redshift bins spanning 
the range from redshift 10 to 3000. 

At each timestep, the photon spectrum is updated with 
the results of the various scattering and pair produc- 
tion processes described in [5], and redshiftcd, and pho- 
tons at sufficiently low energies are tagged as "deposited" 
and removed. The threshold for this "deposition" occurs 
when the photon would on average photoionize an atom 
once per timestep, and as in [5], we choose a timestep of 
c?ln(l + z) = 10~ 3 (it was confirmed in that work that the 
results were converged at such a timestep). Below this 
energy deposition threshold, a prescription for the copi- 
ous secondary electrons from ionization and excitation is 
needed to translate the energy deposition history to an 
ionization history. Such prescriptions have been devel- 
oped in [9~Mllj; for the present problem we restrict our- 
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selves to studying the redshift-dependent partition be- 
tween deposited energy and free-streaming high-energy 
photons. 

In order to improve the speed of the code, we have sim- 
plified the treatment of electron cooling relative to [5]. 
For high-energy electrons (> f keV), inverse Compton 
scattering on CMB photons is by far the dominant cool- 
ing mechanism. As the electron energy falls, the electron 
deposits more of its energy to ionization, excitation and 
heating, but also the energy of photons upscattered by 
the electron falls (as the square of the electron energy), 
and these low-energy photons are promptly absorbed by 
the gas. Thus, in regimes where the initial electron does 
not promptly deposit 100% of its energy, and the parti- 
tion is non-trivial, the vastly dominant energy-loss mech- 
anism is inverse Compton scattering; for the purposes of 
this note, we can therefore ignore the other energy loss 
processes. We have verified that making this approxima- 
tion changes our results at less than the percent level. 

Since we are primarily interested in charge-neutral 
sources of energy injection, we will usually consider in- 
jecting e + e~ pairs rather than electrons alone. This is 
irrelevant for the energy-loss mechanisms important for 
high-energy particles, but when the positrons have cooled 
far enough, they annihilate with ambient electrons, pro- 
ducing gamma rays which can in turn escape rather than 
depositing their energy. We also present results for the 
injection of electrons only; these can be obtained directly 
from the results for e + e _ pairs and photons. Since the 
electron cooling time is many orders of magnitude faster 
than the photon cooling time, it is a good approxima- 
tion that the positron annihilates at the same redshift it 
was injected, so an injected positron is equivalent to an 
injected electron + additional photons at 511 keV and 
below. 

Our results are expressed as a three-dimensional grid 
for each of the particle types, describing the fraction 
of the original particle's energy deposited in any given 
timestep, for a particular redshift-of-injection, and a par- 
ticular initial energy. We will denote the elements of 
this grid as l£*_ fe _ f7 = T e+e - )e -, 7 (4 nj ,4 ep! E k )d\n(l + 
^dcp), where z; n j is the redshift of injection, z^ep is the 
redshift of deposition, and E is the initial energy of the 
particle. 

The files are available online in .fits and .hdf for- 
mat, at http://nebel.rc.fas.harvard.edu/epsilon 

(see Appendix |A| for further details); we also supply a 
Mathematica notebook demonstrating how to read the 
.hdf files and reproduce the calculations in this note. 

III. DISCUSSION 

A. Total energy deposition 

As a first step, let us consider the total fraction of 
the initial particle's energy that is deposited by z = 10, 
^deposited, as a function of initial energy and redshift-of- 



injection. In terms of our grid elements, this is given 
by, 

Pdcpo S itcd(4ij ) E k ) = T e+e- ,e- , 7 ' W 
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Fig. [T] shows the results for electron-positron pairs, pho- 
tons, and electrons unaccompanied by a positron. 

The fraction of deposited energy falls with increasing 
redshift, since the universe becomes more transparent as 
it expands. In general, the fraction of energy eventu- 
ally deposited is higher for particles injected at lower 
energy; this can be attributed to the shorter cooling 
times for such particles, so more of the energy deposi- 
tion takes place at higher redshifts where the universe is 
more opaque. However, electron-positron pairs injected 
very close to threshold deposit a smaller fraction of their 
total energy: their kinetic energy is efficiently deposited, 
but their mass energy is converted into ~ 511 keV pho- 
tons (or equivalently, the mass energy of the positron and 
a single ambient electron is so converted), which at low 
redshifts lie within the "transparency window" and can 
free-stream to the present day. 



B. Effective deposition efficiency and comparison 
with previous results 

Deposition efficiency curves were presented for a range 
of DM models in All of these models represent ex- 
amples of "conventional" DM annihilation, in that their 
annihilation cross section is independent of redshift and 
so the annihilation rate scales simply as the DM number- 
density squared. Suppose that, as in this case, the spec- 
trum of the injected particles is redshift-independent and 
the redshift dependence of the injection rate is known. 
Then integrating our grid over the redshift-of-injection, 
weighted by the appropriate redshift-dependent injection 
rate, yields the total amount of energy deposited at any 
given redshift. 

It is convenient to normalize this energy-deposition his- 
tory by an energy injection history. We will generally 
normalize to the energy injected from annihilation or de- 
cay of the smooth DM component, as appropriate for the 
model we are studying, since this has a simple redshift 
dependence. For example, for DM annihilation and ig- 
noring structure formation, the power injected per unit 
volume at redshift z is given by, 

(^) =(i + *)'nW^ M -, (2) 

\dtdVJ injoctcd m DM 

where (av) describes the annihilation cross section, muM 
is the DM mass, and Odm is (as usual) the cosmological 
DM density as a fraction of the critical density p c . Fol- 
lowing the notation of [5], we can describe the energy- 
deposition history by the dimensionless function f{z), 
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FIG. 1: Fraction of initial energy deposited by 1 + = 10, by a particle or pair injected at a redshift z. The different panels 
show the results for different species. Left panel: electron + positron pair (note "energy" on the x-axis refers to the kinetic 
energy of the electron or positron individually, not the pair - we assume the electrons and positrons have the same spectrum - 
but the plotted values give the fraction of the pair's total (mass + kinetic) energy that is absorbed), center panel: photon, right 
panel: electron (in this case, as there is no way to liberate the electron's mass energy, we plot the fraction of kinetic energy 
which is absorbed). The energy axis on the right panel is truncated because the behavior is uninteresting at lower energies; 
when the electron's kinetic energy is small compared to its mass energy, the deposition fraction « 1. 



defined as, 



dE \ 
dtdV ) 



deposited 
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The f(z) curves describing energy deposition from DM 
annihilation (again ignoring structure formation) can 
readily be extracted from our tables T ijk , for any spec- 
trum of photons and e + e _ pairs produced by the annihi- 
lation. Furthermore, it is straightforward to include cases 
where the spectrum of annihilation products is redshift - 
dependent (for example, in the case where there are mul- 
tiple annihilation channels with different velocity or time 
dependences). 

Defining the interpolating function T(z- m - s , z^ ep , E) as 
above, so T ijk = T(z? nj , z J dcp , E k )d\n(l + z dcp ), we can 
write, 
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wnere dEd\n(i+z) describes the spectrum of injected par- 
ticles (of a particular species), per comoving volume, as a 
function of energy and redshift. I ) refers 

V V ' > J norm 

to the energy injection history we are using for normaliza- 
tion (we will typically pick something simple, like anni- 
hilation or decay of DM neglecting structure formation) . 

For example, studying energy deposition from con- 
ventional DM annihilation, where the energy injection 
histories can be factored into a redshift dependence 
and a spectrum per annihilation dN/dE, it is natural 



(and agrees with the convention in [5]) to normalize to 
the injected energy from DM annihilation in the same 
model, thus canceling out redshift-independent, energy- 
independent proportionality factors in Eq. [4] such as 
(crv). For conventional DM annihilation, 



dN _ f dN \ dt 
dEd\n{l + z) ~ \dEdVdt) d\n{l + z) 

(dN\ {l + zf 
oc - 



\dE J H(z) 



(5) 
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dN cx (1 + z) 6 , dV oc (1 + zf and H{z) 



rfin(i-t-z) , thus yielding, 



H{z) 



species 



^ , {l + z'fdz' f , ^dN 
L / H {z') J T (z,z,E)E — dE. (6) 



species 



For conventional DM decay, with a lifetime much longer 
than the age of the universe (i.e. d Edydt (1 + z ) 3 ) 
and a spectrum-per-decay of dN/dE, the corresponding 
expression is, 



E 



rfln(l + z') 



J T(z / ,z,E)E^dE 



f dN , 

E lV E 



species 



(7) 



Here we have normalized the energy deposition curve to 
the energy injection curve from the same process; in both 
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cases we can therefore think of f(z) as an "effective effi- 
ciency" mapping injection to deposition. 

Fig. [2]shows the efficiency curves for photons and e + e~ 
pairs as a function of injection energy and redshift-of- 
deposition, for injection histories corresponding to anni- 
hilating and decaying DM. In the annihilation case, we 
have included only energy injection corresponding to the 
smooth DM relic density, neglecting the onset of struc- 
ture formation below z ~ 100. 

Note that unlike the total fraction of energy deposited, 
the effective efficiency can rise substantially at low red- 
shifts, particularly in the annihilating-DM case. The in- 
creased transparency of the universe is more than com- 
pensated by the residual pool of photons produced by 
more-rapid energy injection at high redshifts, which re- 
main in the transparency window for a long period and 
only deposit the bulk of their energy at lower redshifts. 
The efficiency is systematically higher for annihilating 
DM than decaying DM, because annihilating DM injects 
more power at higher redshifts when the universe is more 
opaque. 

In Fig. [3] we compare the effective-efficiency curves 
obtained by this procedure to those given in [5] for an- 
nihilating DM, for the same SM final states and DM 
masses listed there. For each final state and DM mass 
we take the injected spectrum of photons and e + e~ 
pairs calculated using Pythia (as in [5 ), integrate these 
spectra over the energy- and species-dependent effective- 
efficiency curves plotted in Fig. [2j and compare to the 
f(z) curves from [5]. We obtain excellent agreement, 
indicating that our rather coarse binning is nonetheless 
adequate for the purpose. 



C. Inclusion of structure formation 

Constraints on DM annihilation from the CMB are 
typically computed only taking into account the smooth 
component of the DM. This is appropriate at high red- 
shift, when density fluctuations are small and perturba- 
tive, and conservative in all cases; it allows a very ro- 
bust and precise constraint, since the cosmological den- 
sity of the DM is constrained at the percent level. How- 
ever, there has been considerable interest in the effects of 
DM annihilation/decay at lower redshifts (z < 30) where 
there might be measurable effects on the 21 cm line of 
hydrogen |12H17j . or where the annihilation byproducts 
might help reionize the universe p"8Tf22] . 

At these redshifts it becomes crucial to account for the 
contribution to annihilation from bound DM halos; since 
the annihilation rate scales as the square of the density, 
annihilation in halos eventually greatly dominates over 
the smooth component. Some previous studies have as- 
sumed an on-the-spot approximation, taking the energy 
deposition rate to be equal to the energy injection rate 
or making some simple ansatz for the fraction of injected 
energy that is promptly absorbed; others have multiplied 
the energy deposition rate (e.g. parametrized by f{z) as 



defined above) by a factor describing the extra annihi- 
lation from halos at that redshift. However, these ap- 
proaches can give rather inaccurate results, because at 
low redshifts the absorption of annihilation products' en- 
ergy is often significantly delayed; the annihilation rate 
in bound halos that is relevant is set by the redshift of 
injection, not deposition. 

There are substantial uncertainties on the contribu- 
tion to the total annihilation rate from bound subhalos 
- the smallest halos are expected to be far below the 
resolution of A-body simulations - but in the present 
context we will use the simple prescription outlined in 
[2"T] to demonstrate the effect of delayed energy depo- 
sition. Employing the Press-Schechter differential mass 
function, the authors of |U] find that the energy injec- 
tion rate from annihilation in halos can be parameterized 
by, 

dVdt \ 1 + z h J m DM 

where the redshift of halo formation Zh and the nor- 
malization factor fh = + zf) 3 /nfw(c/i) should be 
extracted from A-body structure formation simulations. 
[2Tj considers Zh ~ 20 and fh ~ 10 9 — 10 10 to be plausible 
values for these parameters. 

Factoring out the energy injection rate from the 
smooth DM density (Eq. [2]), the enhancement factor 
from including halos is given by, 

1 + -A^erfc ( 1 + 1) . (9) 
(1 + z) 3 \l + z h J 

Given any ansatz for the boost to DM annihilation 
from halos, since we have the full function mapping in- 
jection to deposition, it is trivial to compute the modified 
energy deposition curve: we simply multiply the spec- 
trum of injected particles, dEd $ 1+z) <x 
by the appropriate redshift-dependent halo enhancement, 
inside the integral in Eq. [4| Fig. [4] shows the di- 
mensionless energy deposition function f(z) as a func- 
tion of redshift with the halo parameterization above, 
taking Zh = 20 and fh = 10 9 , normalized to the en- 
ergy injection at that redshift from the smooth compo- 
nent only. (This removes the interpretation of f(z) as 
an effective efficiency, but facilitates comparisons with 
the baseline case where we ignore structure formation.) 
For comparison, in Fig. [5] we show the difference be- 
tween these results and an on-the-spot approximation, 
taking all the injected energy to be promptly deposited 

(i.e. f(z) = 1 + ^j^erfc (j^))- We see that the 
effects of delayed energy deposition can be substantial, 
especially for higher-mass DM. We also show the results 
of the common approximation of simply multiplying the 
f(z) derived for the smooth component by the enhance- 
ment factor from the halos; again, this is a rather poor 
approximation for high-mass DM. 
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FIG. 2: Effective efficiency (i.e. the ratio of the energy deposition history to the energy injection history) for (top) annihilating 
DM and (bottom) decaying DM, for electron-positron pairs (left-hand panels) and photons (right-hand panels), as a function 
of redshift-of-deposition and initial energy of the photon/electron/positron. 



IV. CMB LIMITS ON UNCONVENTIONAL 
INJECTION HISTORIES: A REVIEW 



It is straightforward to obtain deposition histories 
for models with unconventional redshift-dependent en- 
ergy injection histories, simply by changing the factor 
by which we multiply the grid when summing over the 
redshift-of-injcction. Given such a deposition history, we 
can immediately apply the principal-component methods 
of [231 to estimate the CMB constraints on such a class of 



models 1 . For convenience, we briefly review the method 
here. 



Suppose we write the energy deposition history 



dE 
dVdl 



as 



some dimensionless function f(z), multiplying a dimen- 
sionful normalization factor e and a fixed baseline energy 
injection history with no dependence on the model under 



1 These methods include an approximate formulation for the par- 
tition of "deposited" energy between ionization, excitation, heat- 
ing and low-energy (CMB-like) photons for which the universe is 
transparent; these approximations may mis-estimate the bounds 
by several tens of percent 1241 . 
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FIG. 3: Left panel: The f(z) curves derived in [5] for 41 models of interest (solid black lines), and the recovered f(z) curves 
using the results presented in this paper (red diamonds), between z = 10 and z = 1500. Right panel: The percentile difference 
between the two curves, between z = 10 and z = 1500. 
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FIG. 4: The ratio of the energy deposited, from the halos + the smooth component, to the energy injected from the smooth 
component, as a function of injection energy and redshift-of-deposition, for (left panel) e + e~ pairs, (right panel) photons. 



study, b(z). For example, for energy injection from DM 
annihilation with DM mass m x and annihilation cross 
section (crv), we typically choose, 

b(z) = (i + z) 6 n 2 DM c 2 P l 

and for DM decay with lifetime r, 
b(z) = (1 + z) 3 n DM p c c 2 , 



m DM 



(10) 



(11) 



We can now parameterize the redshift dependence of 
the energy deposition, described by f(z), by its coeffi- 
cients in a basis of orthonormal principal components 



{ei} supplied by 



ef{z) 



N 

£ 

i=l 



'ei(z) ■ ei(z) 



ef(z) ■ ei(z), 



where the dot product is defined appropriately as a sum 
over discrete variables or an integral over a product of 
functions. Principal components have been derived by 
[25] for both choices of b(z) described above; of course, 
this is just a choice of basis, but the appropriate choice 
of basis means the effect of the energy injection on the 
CMB can be adequately described by a smaller number 
of principal components. 
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FIG. 5: The ratio of the energy deposited, from the halos 
+ the smooth component, to the energy injected from the 
smooth component, as a function of redshift-of-deposition, 
for three different electron injection energies (assuming they 
are injected as a pair with a positron of the same energy), as 
well as the on-the-spot curve where the deposited energy at 
every redshift is equal to the injected energy (from halos + 
smooth component). Dotted lines show the result of simply 
multiplying the f(z) derived for the smooth component by 
the enhancement factor from the halos, rather than correctly 
inserting the enhancement factor inside the integral in Eq. [4] 



Each principal component &i(z) has an eigenvalue 
Xi associated with it (with different eigenvalues for 
WMAP 7, Planck, a cosmic variance limited (CVL) ex- 
periment, etc), which describes its detectability when 
multiplied by some canonical (dimensionful) normaliza- 
tion factor Si — eVi. The total significance of the sig- 
nal is then given (in the Gaussian approximation) by 
Xief/e, so for example, the 2o limit corresponds 



approximately to the constraint, 



s < 



2e 



Alternatively, e/y/Xi can be replaced by the error bar on 
energy injection corresponding to the ith PC, <7j, as ob- 
tained from a likelihood analysis using CosmoMC, in which 
case the limit is, 



e < 



For Planck forecasting, the two limits are essentially iden- 
tical, as the approximations used to derive the eigenval- 
ues are very good for the energy injections constrained 
by Planck . For energy injections at the WMAP 7 limit, 
one of the approximations begins to break down; more 
accurate results will be obtained by using the CosmoMC 
bound on the first principal component, that is, 



£ < 



1.2 x lfr 26 cm 3 /s/GeV 



Inclusion of ACT data tends to strengthen constraints 
from WMAP 7 by - 15% [25]. 

For Planck and a CVL experiment, where including 
more than one principal component can be justified (see 
[2"3] for details), the values of the first few are, 

• Planck: 

{01,0-2,0-3} = {1.1,2.4,4.1} x lCT 27 cm 3 /s/GeV, 

• CVL: 

{a lt 0-2, 0-3, 0-4, 0-5} = {0.5,1.1,1.8,2.5,3.4} x 
lCT 27 cm 3 /s/GeV. 

Alternative principal components were derived for de- 
caying species (the second b{z) profile given above), since 
decay-like energy injection histories have much greater 
relative effects at low redshift than annihilation-like en- 
ergy injection histories. The principles are identical, but 
a likelihood analysis with CosmoMC has not been per- 
formed. The values of the first few e/\fX~i are: 

• WMAP: 

{0-1,0-2,0-3} = {5.3,6.4,9.5} x 10- 25 s-\ 

• Planck: 

{0-1,02,03} = {1.2,1.6,2.3} x 10- 25 s-\ 

• CVL: 



{0-1,0-2,0-3,04,0-5} 



lO-^s- 1 . 



{2.6,5.7,8.0,12,20} 



V. NEW AND UPDATED CONSTRAINTS 
A. Limits on late-decaying species 

DM itself may decay, with a lifetime much longer than 
the age of the universe; alternatively, new metastable 
states might decay, either to a lighter DM state + Stan- 
dard Model particles (e.g. or purely into Stan- 
dard Model particles (in this latter case, the metastable 
species could only constitute a very small fraction of the 
matter density). If photons, electrons or positrons are 
produced in these decays, they can be constrained by the 
methods discussed here. Such constraints were worked 
out for the 1-year and 3-year WMAP data in [71 H25]. 
employing an on-the-spot approximation. 

We consider a decaying species with lifetime r. Then 
the rate of energy injection per unit volume is given by, 



UtdV J injcctcd 



(1 + z) 3 AQ dccPc c 2 



-t/r 



(12) 



where ^dec is the density of the decaying species (in units 
of the critical density), and A is the fraction of its mass 
energy liberated in each decay. 

As previously, we can normalize the energy deposition 
curve to the profile of energy injection for DM decaying 
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with a long lifetime, b(z) = (1 + z) 3 H-£,mPcC 2 ; specifically, 
we define f(z) by, 

V*^/ deposited V DM T) 

and compute /(z) by summation over the T i]k grid, for 
any spectrum of decay products and lifetime t. (Note 
that now, if deposition precisely traced injection, f(z) 
would be equal to e~*/ T , not 1.) For any value of r, we 
can then constrain Af^dec/^DM by the CMB as described 
above, by taking dot products of f(z) with the principal 
components. The results for electrons and photons in- 
jected at a broad range of energies are shown in Figs. [6] 
and0 

If Afidec = SIdm, i.e. the decaying species consti- 
tutes all the DM, we see that depending on the spec- 
trum of decay products the lifetime is constrained to 
be t > 10 23 - I0 25 s, by WMAP 7 at the 95% confi- 
dence level. For WIMP DM, this constraint is somewhat 
weaker than that obtained from present-day limits [29], 
from which we find r > 10 26 s (albeit with substantial 
dependence on the decay channel). Planck is expected 
to improve the limit from the CMB by nearly an order of 
magnitude. However, the general statement here, that if 
DM is still decaying in the present day then present-day 
probes are more sensitive than the CMB, should not be 
surprising. For t ~ 10 13 — 10 14 s, the limit may be as 
strong as A£7 dcc < 10 _11 fi DM . 

Comparing the forecast constraints for Planck in the 
on-the-spot approximation, our results agree well with 
|28j . except for lifetimes smaller than ~ 10 13 s, where 
our limits are somewhat weaker. This may be due to 
a slightly different treatment of the corrections to re- 
combination: in particular, our constraints are based on 
|23| . which ignored the production of additional Lyman- 
a photons due to excitations from the ionizing parti- 
cles. This is conservative, in the sense that it weakens 
the constraints; the usual alternate treatment (see e.g. 
[Jj) assumes that all energy going into excitations re- 
sults in additional Lyman-a photons, which may over- 
estimate the effect. The difference in the constraints is 
small for conventional DM annihilation (see the discus- 
sion in Appendix A of |23j ) , but this effect is most impor- 
tant around recombination itself: thus we might expect 
it to have the greatest impact in models where the decay 
lifetime is less than the age of the universe at recom- 
bination (since there are no significant constraints from 
decays when the universe was ionized). 



B. Comparison to previous bounds 

There are existing constraints on decaying metastable 
species from big bang nucleosynthesis, large-scale struc- 
ture formation, and present-day indirect detection exper- 
iments. For comparison, we show some of these limits in 
Fig.gj 



Peter et al [31 , 33 35 have studied the case where 
the DM has a nearly-degenerate excited state, which de- 
cays at some late time, providing a velocity "kick" to 
the DM particle. These kicks can modify the evolution 
of DM structure, and so such models can be constrained 
by comparison with observation. There is also a poten- 
tially interesting region where the resulting modifications 
to halos are not ruled out, but do not precisely reproduce 
the results of standard cold DM cosmologies. We show 
both the constrained and the "potentially interesting" re- 
gion, and note that both appear to be ruled out by CMB 
constraints if the decays have a significant branching ra- 
tio to electromagnetically interacting particles. 

Injection of cither leptons or hadrons can modify big 
bang nucleosynthesis, and constrain decaying species 
with lifetimes much smaller than the age of the uni- 
verse at recombination. These constraints depend rather 
strongly on the decay channel; they are much more pow- 
erful for hadronic modes than leptonic modes. We do 
not present an exhaustive survey of such constraints, but 
show one illustrative bound, for a scenario where a mas- 
sive gravitino decays hadronically into neutralino DM 
[3D] . For decay lifetimes comparable to (and longer than) 
the age of the universe at recombination, rather than at 
BBN, the CMB constraints are far more sensitive. 

If the lifetime of the decaying species is comparable to 
the present age of the universe, a wide variety of probes 
from indirect detection experiments become viable; for 
example, studying gamma-rays in the isotropic diffuse 
photon background jMJCH], m the Milky Way [37], and in 
galaxy clusters [29 j . These limits depend rather strongly 
on the spectrum of decay products, and thus on the mass 
of the decaying species and whether all its mass energy 
is liberated in the decay. We show the limits on the 
maximum allowed mass fraction of the DM, as a function 
of lifetime, assuming a particle decaying completely to 
bb or pairs, and scanning over DM masses from 

300 GeV to 10 TeV, using the combined constraints from 
clusters and the extragalactic background [29] . 

Going beyond gamma rays, the limits on decays into 
photons, electrons and neutrinos have been studied by 
[3"2"] . The authors consider the case where the decay liber- 
ates only part of the particle's mass energy; as previously, 
the constraints are highly dependent on the spectrum. 
We plot the most optimistic bound (i.e. choosing the 
mass splitting that leads to the greatest observable sig- 
nal) , scanning over splittings from 4 keV to 60 TeV. This 
approach has advantages over the CMB limits in that it is 
sensitive to decays to neutrinos, but only certain ranges 
of the mass splitting can be constrained (because of the 
sensitivity to the spectrum of decay products). 

In general, we see that if the DM is still decaying in 
the present day, constraints from observation of the de- 
cay products can be much more stringent than the CMB 
bounds (and can even constrain neutrinos, which the 
CMB docs not); however, this statement docs depend 
on the spectrum and species of the decay products. For 
lifetimes smaller than ~ 10 17 s, the CMB constraints 
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FIG. 6: Upper bounds on the fraction of the total DM density that can be composed of a decaying species, as a function of 
the decay lifetime and the energy of the decay products: the 95% confidence limit from WMAP 7 (top row) and the forecast 
95% confidence limit from Planck (bottom, row). The left-hand panels apply to e + e~ pairs (as usual, the energy is for a single 
particle), the right-hand panels to photons. Note that these constraints also apply to the decay of a metastable excited state; 
see the text for details. 



dominate, providing bounds on the maximal mass den- 
sity liberated in decays that are comparable to the best 
present-day constraints. While the on-the-spot approx- 
imation may break down, the constraints vary by less 
than an order of magnitude for photon and electron en- 
ergies ranging from sub-keV to 5 TeV. (Of course, there 
is a model dependence in the branching ratio to these 
electromagnetically-interacting final decay products - the 
limits should be rescaled by this branching fraction - but 
our analysis shows that the dependence on the spectrum 
of the decay products is modest.) 



C. Limits on asymmetric dark matter 



In the "asymmetric dark matter" scenario, the DM 
possesses an asymmetry in the early universe similar to 
the baryonic sector, with the "dark antimatter" annihi- 
lating away completely and the residual DM yielding the 
observed Hdm- These models need not have any DM 
annihilation signature at late times, but one can be gen- 
erated if the symmetric component is repopulated. This 
could occur by decay of a metastable state [35] , or if the 
DM acquires a small Majorana mass which induces oscil- 
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FIG. 7: The 95% confidence limits from WMAP 7 and Planck 
(forecast) on the fraction of the total DM density that can be 
composed of a decaying species, as a function of the decay 
lifetime, for a range of decay products and energies. For each 
lifetime we scan over the constraint for photons and e + e~ 
pairs with (single-particle) energy ranging from 1 keV to 6 
TeV, to obtain the bands displayed; the width of the bands 
reflect the variation in the constraint between different ener- 
gies and decay products. The upper (blue) band corresponds 
to WMAP 7 limits, the lower (yellow) band to Planck. Note 
that these constraints also apply to the decay of a metastable 
excited state; see the text for details. A branching fraction 
less than unity to electrons/positrons/photons will shift these 
bands upward. Values above the lines are excluded. 



lations between the two states [39MTI . In such scenarios, 
annihilation switches back on at some late time; if the 
critical timescale is before reionization, the CMB can ef- 
fectively constrain such models. 

The time dependence of the abundance of the symmet- 
ric component, and hence of the annihilation rate, can be 
quite non-trivial [401 141] ; in scenarios with oscillation, it 
depends on the interplay between annihilation, oscilla- 
tion timescale and scattering of the DM on the baryonic 
matter (which damps the oscillations via decoherence). 
Using the results we have presented here, any arbitrary 
known energy injection history can be constrained, but 
we will not attempt to present a detailed survey of the 
sample annihilation histories in the literature. 

Instead, as an example, we present bounds on a sce- 
nario where the annihilation rate scales as 1 — e~ l l T for 
some timescale t. This time dependence can occur in 
scenarios where the annihilation requires a component 
that is populated by the decay of another species, with 
timescale r; it is also applicable to "oscillating asymmet- 
ric DM" scenarios where DM scattering damps the os- 
cillations, with fractional abundances for the "DM" and 



"anti-DM" states scaling approximately as (1- 



-*/2t 



)/2 



and (1 — e~'/ 2T )/2, respectively. The resulting limits 
may also be approximately applicable to energy injec- 
tion histories with large oscillations in the annihilation 
rate, since sufficiently rapid oscillations in the energy in- 
jection history should cancel out in their effect on the 
energy deposition history, the cosmic ionization history, 
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FIG. 8: As in Fig. [7| but also showing some other con- 
straints from the literature on the parameter space. The blue 
region shows the sensitivity of BBN to a massive gravitino 
decaying hadronically into neutralinos [30]. The purple re- 
gions show the effect of a nearly-degenerate DM excited state 
on structure formation; the dark purple region is ruled out, 
whereas the light purple band has interesting deviations from 
CDM behavior [3T] . The light green band shows the sensitiv- 
ity of the Fermi Gamma-Ray Space Telescope to weak-scale 
DM decaying completely into fi + fjT pairs, using the combined 
constraints from clusters and the extragalactic background 
[29] : all regions above the band are excluded, the width of 
the band corresponds to a scan over DM masses from 300 
GeV to 10 TeV. The light red band shows the same result 
for complete decay into bb pairs. The black lines show the 
best possible current reach of experiments searching for the 
decay products directly, for decay into photons (solid), elec- 
trons (dashed) and neutrinos (dotted), scanning the energy of 
the decay products from 4 keV to 60 TeV [32]. 



and the CMB. However, this should be checked explic- 
itly if precise limits are desired, once the energy injection 
history is known. 

Since we are now dealing with DM particles that can 
only annihilate with their antiparticles, the annihila- 
tion rate per unit volume is given by (<7t;)nDM^DM ~~ * 
(av)n 2 ot / / L in the symmetric limit, where n to t = ^dm + 
n DM is the total DM density. In comparison, for Ma- 



jorana DM the corresponding rate is 



t /2. Since 



the energy injected per annihilation is 2todm, the energy 
injection rate per unit volume is given by, 



dE \ 

injected 

and we define f(z) by. 



— i = (i + zfnl MP 2 c c 2 



(14) 



\dtdV J 



f{z)eb(z), 



deposited 

b(z) 



(i + zfnl MP y, 

2m D M 



(15) 



In the on-the-spot approximation, we would find f(z) = 
1 - e~ t(z ^ T . In the limit where r < Recombination, i-e. 
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FIG. 9: Upper bounds on the parameter for an asymmetric DM model where the annihilation turns on with some 

timescale r, as a function of r and the energy of the decay products: the 95% confidence limit from WMAP 7 (top row) and 
the forecast 95% confidence limit from Planck (bottom row). The left-hand panels apply to e + e~ pairs (as usual, the energy is 
for a single particle), the right-hand panels to photons. 



the DM has oscillated or decayed to a purely symmetric 
phase by the onset of recombination, we should recover 
the results for ordinary annihilating Dirac DM. 



bounds in the small-r limit are, 



Our limits on this scenario are shown in Figs. |9fT0l For 
small r the limit on ^ av ^ asymptotes to a fixed value, as 

todm ... 

expected; for the energy ranges where deposition is most 
efficient (corresponding to / w 1), the 95% confidence 



{< JV ) 
m D M 

m D M 



< 4 x l(T 27 cm 3 /s/GeV (WMAP 7), 

< 9 x l(T 28 cm 3 /s/GcV (Planck). (16) 



These results are in good agreement with those found in 
[23] and g2] (for the WMAP 7 case, they are slightly 
stronger than the limits derived using CosmoMC in [23], 
because the principal component analysis relies on a lin- 
earity assumption which breaks down at the ~ 20% level 
at the WMAP 7 95% confidence limit). Note there is 
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Timescale for annihilation to begin (s) 

FIG. 10: The 95% confidence limits from WMAP 7 and 
Planck (forecast) on the parameter for an asymmet- 

ric DM model where the annihilation turns on with some 
timescale t, for a range of decay products and energies. For 
each timescale we scan over the constraint for photons and 
e + e~ pairs with (single-particle) energy ranging from 1 keV 
to 6 TeV, to obtain the bands displayed; the width of the 
bands reflect the variation in the constraint between differ- 
ent energies and decay products. The upper (blue) band 
corresponds to WMAP 7 limits, the lower (yellow) band 
to Planck. A branching fraction less than unity to elec- 
trons/positrons/photons will shift these bands upward. Val- 
ues above the lines are excluded. 

a factor-of-2 difference in the bound because the con- 
straints in [23] are for Majorana DM, whereas these are 
for Dirac DM; the same value of (av) leads to twice the 
injected power for Majorana DM, and consequently those 
constraints appear to be a factor of 2 stronger. 

VI. CONCLUSIONS 

I have presented the results of a comprehensive calcu- 
lation for the partition between free-streaming photons 
and energy absorbed by the photon-baryon plasma, for 
electrons, positrons and photons injected into the uni- 



verse during the cosmic dark ages. The results apply to 
particles with sub-keV to multi-TeV energies, and can be 
used to immediately compute the energy-absorption his- 
tory and hence the effect on the ionization history and 
the CMB for any arbitrary source of electromagnetically 
interacting particles, with any redshift and energy depen- 
dence. 

These results allow easy computation of constraints on 
oscillating asymmetric DM, decaying species, and mod- 
els for DM structure formation. I have extended previ- 
ous CMB bounds on decaying species by showing how 
the energy-absorption rate varies according to the initial 
energy of injection, and showed that for decays with life- 
times around ~ 10 13 — 10 14 s, Planck will probe models 
where the fraction of the DM mass density liberated in 
decays is as small as 10 -11 . For oscillating asymmetric 
DM, I have taken one sample scenario and shown how the 
constraints weaken as the characteristic timescale for re- 
population for the symmetric component becomes longer 
than the age of the universe at recombination; for early 
oscillations, we recover the usual CMB bounds on DM 
annihilation. When DM structure formation is included, 
the fact that energy is not absorbed at the same red- 
shift at which it is injected becomes very important for 
scenarios with annihilating DM; I have shown how the 
energy of the annihilation products changes the energy- 
absorption history (which will in turn modify the effect 
on reionization). 

The complete set of numerical results used 
in this note, and a Mathematica notebook 
demonstrating their use, are available online (at 
http://nebel.rc.fas.harvard.edu/epsilon) for the 
benefit of the research community. 
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Appendix A: Online Supplemental Material 

We have made the tables 7^+*- and T^ k (see Sec. [n| available at http : //nebel . rc . f as . harvard . edu/epsilon, in 
.hdf and .fits format, labeled as resultsgrid_elec .hdf (.fits) and resultsgrid_phot .hdf (.fits) respectively. 
Each .hdf or .fits file contains six arrays, labeled: 

• OUTPUT _REDSHTFT: this array provides the abscissa for output redshift, i.e. the value of 1 + z at which the energy 
is deposited. 

• L0G10 (ENERGY/EV) : this array provides the abscissa in energy, with values given by log 10 (energy in eV). Note 
that this is kinetic energy; a particle annihilating or decaying to e + e~ would need a mass sufficient to provide 
this energy in addition to the mass energy of the pair. 

• INPUT_REDSHIFT: this array provides the abscissa for input redshift, i.e. the value of 1 + z at which the energy 
is injected. 
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• DEPOSITION_FRACTIONS: this array provides the table T tjk for the appropriate species: that is, for a particle 
injected at some input redshift and energy (given by the abscissa arrays), the fraction of its initial energy 
deposited in the (log-spaced) timestep associated with the output redshift. 

• F_CHECK: this array is redundant and intended only as a cross-check; it describes the effective efficiency /-curve 
for DM annihilation to the species in question, with DM mass given by L0G10(ENERGY/EV) , sampled at the 
redshift points given by the OUTPUT_REDSHIFT array. It can be used to confirm that you are computing the 
effective f-curve correctly from the DEPOSITION_FRACTIONS array. 

• C0NVERSI0N_F ACTOR: this array does not use any of the results of the code, and is defined as H{z)/{\ + z) 3 , 
where 1 + z is given by the INPUT_REDSHIFT array. It is provided for convenience only, as it is the weighting 
factor applied to the DEPOSITION_FRACTIONS array to compute the f(z) curves for DM annihilation (see Eq. |6j). 

For convenience, we also provide .hdf counterparts of the .fits files described in |23) . which are required to 
estimate the CMB constraints (see Sec. [v]), and a Mathematica notebook that reads the .hdf files and performs the 
calculations shown in this article. 
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